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1.  INTRODUCTION 


SIMPLE  EXAMPLES  OP  COMPLEX 
PHENOMENA  IN  PLASTICITY 


by 


Philip  G.  Hodge ,  Jr. 
Professor  of  Mechanics 
University  of  Minnesota 
Minneapolis,  MN  SS4S5 


ABSTRACT 


In  1951  D.  C.  Drucker  presented  a  classic  example  of  a  sim¬ 
ple  three-bar  aystesi  in  which  a  single  monotonical ly- increasing 
external  load  produced  an  actual  plastic  stress  reversal  with 
one  bar  yielding  first  in  tension  and  then  in  compression.  This 
exasiple  is  reviewed  and  built  upon.  Other  simple  examples 
Illustrate  general  loading,  unloading,  and  reloading  of  elastic/ 
plastic  structures,  shakedown,  and  lack  of  uniqueness  in  contained 
plastic  deformation. 
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There  are  many  complications  in  applications  of  plasticity 
theory,  particularly  in  comparison  with  the  linear  theory  of 
elasticity.  If  thesa  complications  are  encountered  for  the  first 
time  in  s  "real"  problem  with  large  amounts  of  analysis  and  com¬ 
putation  even  for  an  elastic  material,  the  peculiar  features  of 
plastic  behavior  may  be  masked,  and  the  probability  of  undetected 
errors  is  greatly  increased,  ror  this  reason,  there  is  an  obvious 
advantage  to  discussing  plasticity  first  in  terms  of  simple 
structures,  the  simpler  the  better.  Of  course,  if  s  problem  is 
too  simple,  then  some  essential  features  of  plasticity  may  not 
appear.  However,  one  of  the  surprising  features  of  plasticity 
theory  is  that  so  many  of  its  important  features  can  be  Illustrated 
with  problams  which  are  easily  comprehended  by  students  in  a 
beginning  class  In  statics. 

In  the  present  paper  we  shall  illustrate  several  of  the 
features  of  plasticity  theory  with  two  very  simple  three-bar 
statistically  indeterminate  trusses,  as  shorn  in  Figs.  1  and  8. 
Specifically,  in  Sec.  2  we  will  use  the  truss  of  Pig.  1  to  show 
the  difference  between  hardening  and  perfectly-plastic  materials, 
ana  the  effect  of  including  or  neglecting  elastic  strain  components. 
Section  1  extends  thssa  results  to  unloading  and  reloading.  He  then 
examine  some  undesirable  effects  which  may  occur  when  loads  are 
applied  repeatedly.  The  next  two  sections  illustrate  types  of 
behavior  which  are  drastically  different  from  those  encountered  in 
elasticity.  In  Sec.  S  we  look  at  ths  truss  in  Fig.  8  and  show  that 
under  s  monotonically  increasing  single  load  s  bar  may  yield  first 
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in  tension  and  then  In  compression .  Thus  proportional  • loading" 
does  not  by  any  means  guarantee  proportional  "stressing".  Section 
6  returns  to  the  truss  in  Fig.  1  and  shows  that  a  boundary-value 
problem  that  is  "well-posed"  in  the  theory  of  elasticity  may  have 
a  non-unique  solution  in  plasticity.  The  final  section  of  the 
paper  will  quote  analytic  and  numerical  solutions  to  more  complex 
problems  where  the  same  types  of  responses  are  observed. 

2.  LOADING  OP  DIFFERENT  PLASTICITY  MODELS 

The  defining  equations  for  any  problem  in  continuum  or  struc¬ 
tural  mechanics  come  from  three  sources t  statics,  kinematics, 
and  constitutive.  For  the  truss  in  Fig.  I,  the  two  statics 
equations  are  obtained  from  horlsontal  and  vertical  equilibrium  of 
point  D: 

(Tj  -  r3i//I  -  H  (Tj  ♦  r3i//I  ♦  Tj  -  »  (1) 

where  F A  is  the  tensile  force  in  bar  l.  The  kinematics  expresses 
the  elongation  e^  of  each  bar  In  terms  of  the  displacements  u  and 
v  of  point  Ds 

ej  -  (u  ♦  vl//7  e2  -  v  e3  -  l-u  ♦  v)//I  (2) 

Obviously,  we  are  using  a  strictly  linear  theory  in  writing  Eqs. 

(1)  and  (2)  . 

in  the  preeent  section  we  ere  considering  only  the  ceee  where 
stress  end  strain  srs  both  increasing  functions  of  time,  ao  that 
we  can  quite  generally  writs 

o  -  f(c)  (J) 

Therefore,  for  each  bar  of  the  truss  we  may  write 


Fi  *  *lfiUl,/Li  Ml 

where  A^  is  the  cross-sectional  area  and  Lj  la  the  length  of  bar  i. 
Although  f(c)  could  represent  an  experimentally-obtained  stress- 
strain  curve,  we  will  consider  here  the  four  piecewise-1 inear 
curves  model  shown  in  Fig.  2.  Figure  2a  portrays  an  "elastlc/atrain- 
hardening"  material  (E/SH)  whose  equations  may  be  written 


0  <  o  <  oy  a  «  Et 

oy  <  o  o  «  g't  ♦  |1  -  EVE)oy  ($) 

where  E  is  the  elastic  modulus,  P*  is  the  plastic  modulus,  and  oy 
la  the  yield  stress.  Me  define  the  yield  force,  elastic  stiffness, 
and  plastic  stiffness  of  bar  i  by 


"iV 


respectively,  to  obtain  the  constitutive  equation  for  bar  1* 


K/  0  <  Fi  *  Yi  ri  "  ki*i  (7a) 

B/8H  Yt  <  Ft  rA  -  ki*e1  ♦  Y^l  -  kj'/hj)  (7b) 

Similarly,  the  constitutive  equation  for  the  "elastic/perfectly- 
plastlc"  materiel  (B/PP)  in  Pig.  2b  conalste  of  Eq.  (7a)  for 
0  <  rt  <  Yt  end 

B/PP  eA  >  Yj/kj  PA  -  ya  <7c) 

For  the  "rigld/strein-hsrdening"  (R/SH)  and  “r igid-per fectly-plaatic* 
(R/PP)  meteriale  in  Figs.  2c  and  2d,  respectively,  the  equations  are 


R/ 

R/SH 

R/PP 


-  0 


Fi  *  Vi  *  ¥i 
Fi  *  Yi 


(7d) 

(7e) 

(7f) 
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At  this  point  the  per fectly-plastlc  trues  cannot  tolerate  any 
further  Increase  in  load,  but  can  continue  to  deform  indefinitely 
with  no  further  load  Increase,  Thus  the  aolution  for  Stage  3 
consists  of  Eqs.  (17a)  with  any  value  of  v  greater  than  2Y/k. 

The  dashed  curve  in  Fig.  1  shows  the  load -def onset ion  curve  for  an 
E/PP  material  as  defined  by  Eqs.  (11,  12.  It,  17). 

If  the  truss  is  made  of  a K/8H  material,  it  will  not  be  able 
to  deform  at  all  until  all  three  bars  hava  reached  yield.  Thus 
Stages  1  and  2  of  the  previous  models  are  condensed  into  s  single 
Stage  1  with  •  The  bar  forces  are  non-unique  until  Stage  IL 

is  reached  when  they  are  given  by  Eqs.  (17s)  with  v  -  0.  For  any 
positive  v,  then,  all  three  bars  sra  plastic,  and  the  solution  la 
Stage  2,  all  plastic  -  fj  -  Y  ♦  k'v/2  Fj  •  Y  ♦  k*v 

v  -  k'v(i//i  ♦  d  ♦  i/i  ♦  m  i 

The  dot -dashed  curve  in  Fig.  3  shows  the  resulting  load-deflection 
curve. 

Finally,  the  E/PP  model  predicts  sero  displacement  until 
Stage  IL  as  given  by  Eqs.  (17e)  with  v  -  0.  For  Stage  2  the  static 
quantities  are  still  9lven  by  Ills)  with  any  positive  v.  The  load- 
deflection  curve  is  shown  dotted  in  rig.  3. 

Figure  3  is  drawn  with  dimensionless  scales  with  an  assumed 
ratio  k/k’  -  0.1.  Let  us  define  thg  "yield-point  load*  of  the 
truss  as  that  value  V  -  vQ  for  which  the  truss  becomes  very  much 
more  flexible.  For  the  E/SH  model  it  is  given  by  the  last  Eq.  (14), 
and  for  the  other  three  models  it  has  the  precise  value  ¥(/7  ♦  1) . 
Therefore,  if  only  v0  is  required.,  even  the  simplest  R/PP  modal 
gives  a  very  good  approximation  to  It. 


For  values  of  the  load  less  tna*  the  yield-point  load,  the 
effect  of  strain  hardening  is  very  small,  and  it  appears  quite 
reasonable  to  use  the  B/PP  model  as  an  approximation  of  the  more 
accurate  E/SH  one.  On  the  other  hand,  if  the  applied  load  is  to 
be  much  larger  than  V0,  then  the  effect  of  elastic  strains  becomes 
increasingly  less  important,  and  the  r/SH  model  gives  a  good 
approximation. 

3.  UNLOADING  AND  RELOADING 

Equations  (5)  are  no  longer  adequate  to  describe  the  con¬ 
stitutive  behavior  If  the  stress  in  s  bar  is  allowed  to  both 
increase  and  decrease.  Instead,  we  must  formulate  them  in  terms 
of  rates  or  Increments.  To  this  end,  at  *ny  given  instant  we 
characterise  the  behavior  of  a  bar  as  either  elastic  (E) .  strain- 
hardening  (SH) ,  or  parfectiy-plastic  (PP)  .  Then  the  rate  equations 
corresponding  to  (S|  may  be  written 


Bi 

0 

•  Ec 

(19a) 

8Kt 

0 

•  E*  c 

(19b) 

PPi 

o 

-  0 

(190 

Equations  (II)  do  not  form  a  complete  description  of  material 
behavior,  because  we  still  require  rules  to  tell  which  set  ot 
equations  are  applicable  at  any  given  time.  For  the  PP  material 
this  description  Is  quite  simple.  The  magnitude  of  the  stress  can 
never  exceed  the  initial  yield  stress  If  o  -  oy  and  t  is 

increasing  or  if  o  •  -o^  and  t  is  decreasing,  then  the  bar  is 
plastici  otherwise  it  is  elastic. 
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where  we  have  defined 

r  •  k*/k  (23) 

In  stage  4  we  decrease  v  which  will  cause  all  bars  to  revert 
elastic  behavior  until  bar  2  yields  in  compression.  Por  KH  the 
total  change  in  F2  during  stage  4  will  be  AFj  *  -2Y.  Therefore, 
using  (20a)  for  all  three  bars  we  may  write  the  total  change  In 
solution  during  stage  4  as 

A r/Y  -  Avk/Y  -  -2  APj/Y  -  1  (24) 

Adding  these  increments  to  the  velues  at  Stage  3L,  we  obtain 
Stage  4L(RH)i  kv/Y  -  6  Pj/Y  -  3r  Pj/Y  -  -1  ♦  7r  (25) 

rn  stage  5  bar  2  la  plastic  so  we  use  (20a)  for  bar  1  and 
(20b)  for  bar  2.  This  stage  will  and  when  bar  1  becomes  plastic, 
so  the  total  increments  are 

APj/Y  -  -1  Avk/V  -  -2  APj/Y  •  -2r  (26) 

Addition  of  Eqs.  (25)  and  (26)  then  gives 

Stage  5L0CH) i  kv/Y  -  4  Pj/y  -  -1  ♦  Jr  Pj/Y  -  -1  ♦  5r  (27) 
Por  the  rest  of  the  unloading  both  bars  are  plastic  and  take 
Eq.  \ 20b) .  This  stage  ends  when  we  again  revets#  v  and  start  to 
reload.  During  reloading  we  enter  stages  7,  I,  and  9  corresponding 
to  sli  bars  elastic,  bar  2  plastic-  and  all  bars  plastic  respectively. 
The  various  solutions  era 


Stage  6LOCH)  t 

kv/Y  -  -6 

P1/Y  «  -1  -  2r 

F2/Y  -  -1  -  5r 

(2Ba) 

Stage  7’fKH)] 

kv/Y  -  -4 

Pj/Y  -  -2r 

P2/Y  -  1  -  Sr 

(29b) 

Pj/Y  -  1  -  2r  P2/Y  -  1  -  3r 
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Stage  9L(KH) >  kv/Y  -  0  Pj/Y  -  I  -  r  Pj/Y  -  1  -  r  (28d) 

The  load  in  each  stage  is  given  by  the  second  Eq.  (10) . 

Por  an  1H  material  the  aolutlon  will  be  the  same  up  through 
Stage  3L,  and  the  truss  will  have  the  same  rates  in  Stage  4. 

However,  we  must  now  keep  track  of  the  current  yield  force  Y*  for 
each  bar.  During  Stage  2  '  will  increase  and  remain  equal  to  P^, 

and  during  Stage  3  both  yield  forces  will  increase  with  the  bar 
forces.  Thus  at  Stags  3L  the  yield  forces  will  be 

Stage  3L(IH) t  Yj'  •  1  ♦  3r  Yj'  •  1  ♦  7r  (29) 

and  these  values  will  continue  to  hold  in  Stage  4.  Stage  4L  will 
occur  when  Pj  -  -Y^* ,  hence 

Stage  4MIH)i  Pj/Y  4r  Fj/Y  -  -(1  ♦  7r)  kv/Y  -  6  -  l«r  (10) 

During  Stage  5  Y^*  will  Increase  with  the  magnitude  of  ,  but  Yj* 
will  still  be  given  by  the  first  Eq.  (30).  At  the  end  of  Stage  5 

we  have 

Stage  5L( IH) i  F./Y  -  -|1  ♦  3r)  kv/Y  *  4  -  I2r 

1  -  (JD 

F2A  •  -U  ♦  9r  -  2r  ) 

Continuing  in  this  fashion  we  obtain  the  solutions  for  the  rest 
of  the  prescribed  program.  The  results  are  summarised  in  Table  1. 

If  r  >  0.0674  for  an  IN  matarial,  the  yield  force  in  bar  l  will 
incraaea  to  tha  point  that  v  -  0  during  raloading  while  bar  1  is 
still  alsstlc. 

The  results  for  s  PP  materiel  can  be  obtained  from  either  of 
the  above  solutions  by  setting  r  •  0,  or  can  easily  be  found 
directly.  The  results  are  summer! red  in  Table  2. 


Stage  BLOCH)  t  kv/Y  -  -2 


( 2Bc) 


Por  all  solution*  the  external  load  V  is  found  by  substituting 
th*  listed  bar  forces  in  ths  second  Eq.  (10).  rigure  5  shows  the 
resulting  load -displacement  curves  for  the  three  different  materials. 
Results  for  the  rigid/plastic  materials  will,  of  course,  be  very 
similar. 

4 .  SHAKEDOWN 

We  consider  first  an  E/PP  material  which  is  subjected  to  a 
prescribed  loading  program  which  continues  indefinitely,  usually 
as  a  repeated  cycle.  It  the  iomda  are  alw ays  sufficiently  small, 
the  structure  will  remain  everywhere  elastic,  and  Is  of  no  concern 
in  th*  present  context.  At  th*  other  extreme,  if  th*  loads  ever 
exceed  the  yield-point  load  of  the  structure,  there  will  be  no 
equilibrium  solution.  The  structure  will  collapse  catastrophically, 
and  the  rest  of  the  loading  program  will  be  meaningless.  In  the 
present  section  w*  are  concerned  with  loading  programs  between 
these  two  extremes,  i.e.,  the  loading  program  includes  loads  which 
exceed  the  elastic  limit  but  which  nowhere  exceed  the  yield-point 
load . 

W*  illustrate  various  loading  programs  with  regard  to  the 
truss  in  Pig.  1.  in  every  case  we  take  all  three  bars  to  have  th* 
same  cross-sectional  area  and  th*  same  materiel  properties. 

Consider  first  a  program  in  which  H  -  0,  always,  and  V 
oscillates  between  0  end  2.2Y.  During  th*  initial  increase  of  V 
Bqa.  (21).  (12),  and  (It)  apply  but  stage  2  ends  with  the  limiting 
value  of  V* 


Stage  2Li  Fj  -  P3  -  0.6/7Y  P2  -  Y  kv/Y-  1.2/7  V-2.2Y  (12) 

When  V  is  decreased  from  2.2Y,  the  changes  in  both  bars  will  be 
elastic.  Therefore,  it  follows  from  (20a)  and  the  second  Eq.  (10) 
that 

APl-0.5kftv  AFj-kAv  AV  -  ( 1  ♦  1//I)  kAv  (3)) 

Setting  AV  •  -2.2Y  o  reduce  V  to  tero,  leads  to  the  values 
Stage  3l:  Fj  -  Pj  -  ( 1 . 7/1  -  2 . 2) Y  F^  ■  (2.2/7  -  3. 4)  Y 

(34) 

kv/Y  -  3.4/7  -  4.4 

Since  both  forces  are  within  the  elastic  range,  this  solution  is 
the  valid  one  for  stage  3L.  Purther,  when  v  is  again  increased 
to  2.2Y,  Eqw.  (33)  will  hold  for  th*  entire  reloading  process, 
and  the  solution  at  stage  4L  will  be  exactly  the  ism  as  that  at 
Stag*  2L  as  given  by  Eqa.  (32).  Clearly  further  cycles  will  simply 
alternate  between  Eqa.  (32)  and  (34)  with  both  bars  continuing 
to  behave  elastically. 

The  above  behavior  can  be  summarised  as  follows.  There  is  a 
limited  amount  of  plastic  flow  at  the  beginning  of  the  loading 
process,  but  after  this  has  taken  place  the  rest  of  the  cycle  is 
completely  elastic.  A  process  with  this  property  is  called  a 
"shakedown”  process. 

As  an  example  of  a  process  which  doss  not  shake  down,  suppose 
thet  after  the  Initial  loading  to  2.2Y  the  load  v  is  alternated 
between  -2.2Y  and  +2.2Y.  Stages  1  and  2  will  be  the  same  as 
before,  and  the  change  in  Stage  3  will  still  be  given  by  Eqs.  (33) . 
However  this  solution  is  valid  only  until  A P^  *  -2Y,  which  leads  to 
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Stage  3L;  F.  •  (0.6/I-DY  Fj  -  -Y  kv/Y«l. 2/1-2  V--(0.4*/5)Y  (35) 

Further  unloading  takes  place  with  bar  2  plastic,  hence  A F^  -  0. 
Choosing  AV  to  bring  V  to  its  final  value  of  -2.2Y,  using  Eq.  (20a) 
for  AFj ,  and  the  second  (10)  for  V,  we  obtain 

Stage  4L:  F1»-0.3/Iy  Fj  -  -Y  kv/Y  - -0.6/1  V--2.2Y  (36) 

The  reloading  process  is  similar,  consisting  of  Stage  5  with 
all  bars  elastic  and  Stage  6  with  bar  2  plastic  in  tension.  The 
results  are  summarized  in  Table  3a.  We  observe  that  the  solution 
for  Stage  6L  is  exactly  the  same  as  the  one  at  Stage  2L,  so  that 
this  cycle  of  4  stages  will  be  exactly  repeated  in  each  cycle  of 
V  from  +2.2Y  to  -2.2Y  and  back  to  +  2.2Y  as  given  in  Table  3a. 

In  particular,  the  value  of  the  displacement  will  never  exceed 
1 . 697Y/k  which  is  less  than  double  the  maximum  elastic  displacesMnt. 

At  first  glance  it  might  appear  that  this  loading  program 
was  a  safe  one  for  the  truss  since  bar  1  never  yields  and  the 
displacement  is  always  bounded.  However,  from  a  materials  point 
of  view,  the  process  of  plastic  dlsplacesient  is  not  e  reversible 
one.  Although  the  gross  displacement  of  bar  2  le  the  same  after 
a  complete  cycle,  the  tensile  plastic  behavior  in  Stage  6 
does  not  "undo*  the  compressive  plastic  elongation  in  Stage  4, 
but  rather  superimposes  a  tensile  plastic  elongation  on  it.  To 
obtain  some  insight  as  to  why  this  behavior  is  undesirable,  we 
introduce  the  concept  of  "plastic  work". 

The  differential  internal  work  done  on  a  bar  la  given  quite 
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elongation,  so  that  all  work  done  is  plastic.  Further,  since  the 
force  has  its  constant  yield  value,  and  since  the  elongation  change 
is  positive  in  tension  and  negative  in  compression,  we  can  write 
the  Increment  of  plastic  work  in  the  form 

4»p  -  Y  | A, |  (371 

In  particular,  since  bar  2  is  the  only  one  to  yield  in  the  present 
example,  we  may  use  the  middle  Eq.  (2)  to  write 

IY(Av|  when  bar  2  plastic 

(38) 

0  when  bar  2  elastic 

Equation  138)  shows  that  plastic  work  is  done  during  each  even- 
numbered  stage  and  la  always  positive.  The  last  column  of  Table  3a 
shows  the  cteaulative  plastic  work  through  Stage  6.  Clearly  a 
further  increment  of  1.394  Y2/k  must  be  added  during  each  additional 
half  cycle.  Since  the  capacity  of  any  real  material  to  absorb 
plastic  work  is  limited,  the  truss  will  become  unserviceable  after 
e  relatively  email  number  of  cycles.  The  solid  curve  in  Fig.  6 
shows  the  accumulation  of  plastic  work  during  the  first  few  cycles. 

A  loading  program  in  which  this  phenoaten  occurs  is  called  "alternating 
collapse". 

In  Sec.  2  we  concluded  that  it  was  reaeonable  to  neglect 
hardening  at  loads  less  than  the  yield-point  load,  at  least  for 
a  single  loading.  To  what  extent  le  that  conclusion  valid  in 
relation  to  shakedown?  It  turne  out  that  the  answer  is  quite 
different  depending  upon  the  type  of  hardening.  M«  consider 
first  an  B/RH  materiel.  As  before,  bar  1  is  always  elastic,  so 


generally  by  -  F  de.  For  an  B/PP  bar  at  yield  there  ie  no 
change  in  the  bar  force  and  hence  no  change  in  the  elastic 
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we  will  be  dealing  with  only  two  different  incremental  solutions: 


APj  -  kAv/2 

always 

(39a) 

AF2  ■  kAv 

AV  -  kAv ( 1  ♦  1//?) 

bar  2  elastic 

(39b) 

AP2  -  k  *  Av 

V  -  kAv(r  ♦  1//2) 

bar  2  plastic 

(39c) 

During  the  initial  loading  (39b)  will  apply  with  »  f  to 
determine  Stage  1L,  then  (19c)  should  be  used  with  AV  -  (1-2  -  0.5/2)* 
for  Stage  2L.  Stage  3L  is  determined  froei  (39b)  with  AFj  «  -2Y, 
and  Stage  4L  from  (39c)  with  AV  *  (-2.4  ♦  /7)V.  Similar  reasoning 
applies  toStagesS  and  6.  The  results  are  summarised  in  Table  3b. 
Since  conditions  are  exactly  the  same  at  Stages  6L  and  2L,  further 
cycles  will  simply  repeat  the  last  four  lines  of  Table  3b. 

When  a  bar  is  plastic,  the  incremental  work  can  be  written 

du-Fda-Fadt-  Fl/k'  dt  -  d(FJ/Jk‘l  (60) 

However,  to  find  the  plastic  work  differential  we  suit  subtract  the 

elastic  differential  dw  -  d(P^/2k) .  If  and  AF  represent  the 

beginning  value  and  increment,  respectively,  of  F  during  a  plastic 

stage,  we  may  write  the  plastic  work  Increment  in  the  form 

4N  -  O.SU/k1  -  l/k||(F„  ♦  4F|J  -  F„J| 

P  Q  U  (4U 

-  (1  -  rl (4F/k*»  (,n  .  JF/2) 

In  particular,  for  bar  2  Bq.  (41)  becomes 

AWp-(l  -  r)Av(F0  ♦  0.5AF)  (42) 

The  last  two  columns  of  Table  3b  show  the  increments  and  cumulative 
values  of  the  plastic  work,  and  the  dashed  curve  in  Fig.  6  shows 
the  cumulative  plastic  work  over  the  first  few  cycles  for  r  -  0.1. 

Comparing  the  results  for  the  8/PP  and  E/KH  models  we  see  that 
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the  former  is  a  reasonable  approximation  to  the  latt  parti¬ 

cular,  both  ax»dels  predict  that  the  plastic  work  wil 
without  bound  as  the  cycles  continue. 

The  B/IH  material  behaves  quite  differently.  The  increments 
are  still  governed  by  Bqs.  (39),  but  the  passive  yield  force  comes 
from  (19h)  rather  than  from  ( 19 i )  .  The  resulting  solution  is 
exactly  the  same  through  Stage  3,  but  Stage  3L  occurs  when  j» 
equal  to  the  negative  of  its  value  at  Stage  2L,  rather  than  when 
AF  -  -2*.  The  net  result  is  that  the  truss  spends  more  of  its 
unloading  time  in  the  elastic  Stage  3  and  less  in  the  partxally- 
plastic  Stage  4.  The  Some  thing  happens  during  the  reloading 
phase.  Therefore,  ss  shown  in  Table  3c,  the  solution  at  Stage  feL 
is  different  from  that  st  Stage  2L,  so  that  the  next  cycle  of 
unloading-reloading  will  be  different.  In  particular,  we  note 
that  the  value  of  F^  at  the  end  of  each  cycle  has  increased. 

Since  this  value  is  the  current  yield  stress  of  bar  2,  a  larger 
proportion  of  each  succeeding  cycle  is  spent  in  s  fully  elastic 
state. 

Results  for  the  first  few  cycles  are  shown  in  Table  3c  and 
by  the  dotted  curve  in  Fig.  6.  Although  bar  2  is  plastic  during 
part  of  every  cycle,  it  is  apparent  that  the  process  is  con¬ 
verging  to  an  entiraly  elastic  ona ,  and  that  tha  total  amount  of 
plastic  work  is  finite. 

Let  us  return  to  the  E/PP  material  and  consider  a  different 
loading  cycle.  We  begin  In  the  same  way  by  increasing  V  to  2.2Y, 
but  we  then  superpose  en  alternating  load  at  45*.  Therefore, 
Stages  1  and  2  are  the  same  as  in  the  previous  example,  but 
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beginning  with  Stage  3  the  load  increments  are  defined  by 
AH  =  -AV  -  AX//2,  where  X  is  cycled  between  0  and  1.1Y. 

Therefore,  F.gs .  (1)  can  be  written 

f.Pj  -  AFj  -  4»  APj  *  4Fj  +  /lAFj  *  -4»  (43) 

At  the  end  of  Stage  2  bar  2  is  plastic,  but  the  upward  component 
of  the  diagonal  load  reduces  P2  so  that  Stage  3  is  fully  elastic. 

The  increments  of  the  three  bar  forces  and  two  displacement  com¬ 
ponents  are  all  proportional  to  AX  and  are  given  in  the  first 
column  of  Table  8  .  Although  bars  1  and  3  are  both  in  tension  at 
the  start  of  Stage  1,  Fj  is  increasing  and  F^  is  decreasing 
durinq  the  stage.  Stage  It  is  determined  by  bar  1  reaching  its 
tensile  yield  force  when  *  •  (0.8  -  0.2/7JY  «  0.517Y. 

In  Stage  4  AFj  *  0  and  Egs.  (7a)  are  valid  for  bars  2  and  3, 

Toq ether  with  (83)  this  leads  to  AF^  -  0,  also.  The  complete 

solution  is  shown  in  column  2  of  Table  8.  Stage  4  terminates 

when  x  equals  its  final  value  of  1.1Y.  The  complete  solution  for  i 

Stage  8L  is  given  in  Table  5. 

Me  now  reduce  the  diagonal  load  so  that  A A  is  negative.  This  j 

causes  a  decrease  in  Fj ,  so  that  the  truss  is  again  fully  slsstic. 

As  shown  in  Table  8,  a  negative  AX  increaaes  F^  and  bar  2  will 
yield  in  tension  when  A  ■  0.583.  The  complete  solution  for  the 
resulting  Stage  5L  is  shown  in  Table  5.  In  stage  6  AF^  “  0 
which  leads  to  the  same  force  increment  values  as  In  Stage  4 
(see  Table  4).  However,  bars  1  and  3  ara  now  elastic  and  the 
corresponding  Rqs.  (7a)  show  that  Au  ♦  Av  -  0,  whereas  in  Stage  4 
Av  -  0.  Thus,  the  displacement  increments  ara  different  in  the 
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two  stages.  Stage  6L  ia  defined  by  1  -  0,  and  the  resulting 
complete  solution  is  given  in  Table  5. 

The  bar  forces  are  all  exactly  the  same  at  Stage  6L  as  they 
were  at  Stage  2L.  Therefore,  they  will  repeat  the  same  sequence 
of  values  during  each  additional  cycle  of  load.  However,  the 
vertical  and  horizontal  displacements  have  each  increased  by 
(0.3  ♦  0.2/2)Y/k  »  0.583Y/k  during  the  cycle.  Clearly  they  will 
increase  by  the  same  amount  during  each  succeeding  cycle.  There¬ 
fore,  although  at  any  given  instant  the  truss  has  at  least  two 
bars  elastic  and  can  support  the  given  load,  the  deformations 
grow  indefinitely  with  time  and  the  truss  eventually  becomes 
unserviceable.  This  behavior  has  been  termed  "»nci  mental 
collapse".  A  history  ->f  the  displacement?  during  the  first  few 
cycles  is  shown  in  Fig.  7. 

Physically,  we  can  observe  that  if  bars  1  and  2  were  both 
in  plastic  tension  at  the  same  time  the  truas  would  be  a  mechanism 
capable  of  a  rotation  about  the  end  of  bar  3  with  u  •  v.  Of 
course  this  does  not  occur,  but  in  the  course  of  a  complete  cycle 
bars  1  and  2  are  plastic  at  different  times.  The  net  effect  at 
the  end  of  a  complete  cycle  is  a  limited  mechanism  motion 
Au  ■  Av  »  0.583  Y/k,  as  shown  by  the  detailed  analysis. 

For  this  loading  program  the  two  stainhardening  models  are 
the  same.  Stage  2*,  is  again  given  by  the  appropriate  line  from 
Table  3b  or  3c.  From  then  on  the  equilibrium  equations  are 
given  by  (43)  .  Combining  the  kinematic  end  constitutive  equations 
we  obtain 


I 
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E/  AFj  “  kUv  ♦  Au)/2  AF2  -  kAv  AF3  -  k(Av  -  Au)/2  (44a) 

/SH  AFj  *  k* (Av  ♦  Au)/2  APj  ■  k'Av  AFj  -  k 1 ( Av  -  Au) /2  (44b) 

Stage  3  is  fully  elastic  and  terminate*  when  bar  1  reaches  its 
yield  force,  Y.  The  complete  solution  for  this  and  subsequent 
atages  are  given  in  Table  6.  In  Stage  4  we  use  (44b)  for  bar  1 
and  (44a)  for  bars  2  and  3,  and  terminate  the  stage  when  X  reaches 
its  max  mu*  value,  1.1Y. 

Decreasing  X  in  Stage  5  causes  all  bars  to  becoaie  elastic  again. 
However,  the  stage  does  not  end  until  bar  2  reaches  its  current 
yield  stress,  V^' ,  This  value  is  the  saxisua  of  previous  values 
of  F2,  i.e.,  Y2*  »  1.061  which  was  attained  at  Stage  2L. 

Stage  6  continues  with  bar  2  plastic  and  ends  when  X  is  reduced 
to  0.  Since  Stages  6L  and  2l  do  not  have  the  same  bar  forces,  if 
we  want  More  stages  we  Bust  find  the*  explicitly.  The  qualitative 
description  of  each  cycle  is  the  saw:  as  X  is  Increased  all  bars 
are  elastic  at  the  start  but  bar  1  become*  plastic  before  X  reaches 
1.1Y;  as  X  is  decreased  all  bars  are  first  elastic  but  bar  2 
becomes  plastic  before  the  cycle  ends.  However  each  time  bar  1  or 
bar  2  yields  its  yield  force  is  Increased  so  that  during  the  next 
cycle  the  truss  will  regain  fully  elastic  until  nearer  the  end  of 
the  cycle.  Table  6  gives  the  complete  results  through  Stage  14L, 
corresponding  to  three  complete  cycles. 

Of  particular  interest  are  the  total  displaceaients  added 
during  each  cycle,  using  the  results  of  Table  6  it  is  easy  to 
show  that  the  ratio  of  the  increase  in  either  u  or  v  to  its 
increase  in  the  previous  cycle  is  the  constant  value  n  “  0.571636. 
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Thus  the  net  displaceawnt  can  be  written  as  a  geometric  series 
which  converges  to  the  final  values 

u  -  1 . 508Y/k  V  -  2.932Y/k  (45) 

Therefore,  at  least  in  this  simple  example,  the  introduction  of 
strainhardening  leads  to  finite  values  of  the  total  displacements 
as  compared  with  the  Infinite  values  predicted  by  the  E/PP  model. 

5.  REVERSE  STRESSING  UNDER  HO NOTONIC  LOADING 
In  the  previous  sections  we  have  seen  various  examples 
where  the  force  ratios  between  bars  will  change  as  one  or  mote 
bars  becomes  plastic,  even  though  the  external  load  ta  uniformly 
increasing.  Here  we  present  a  simple  example  in  which  the  tatios 
not  only  change  magnitude,  but  actually  change  sign.  The  truss 
in  Fig.  Ba  consists  of  three  numbered  vertical  bars  made  of  an 
E/PP  material  and  joined  to  a  perfectly  rigid  horizontal  bar. 

The  three  vertical  here  have  equal  areas  A,  lengths  L  and  moduli  E, 
but  are  assumed  to  have  different  yield  stresses  given  by 

Yj  -  Y  Y2  -  20Y  Yj  -  1 2Y  (4b) 

A  vertical  load  p  ia  applied  as  indicated.  It  is  increased 
just  to  its  yield-point  value,  and  then  decreased  to  zero.  There 
are  no  horizontal  loads,  and  we  assume  zero  horizontal  motion. 

Static  equations  are  obtained  by  considering  vertical  and 
moment  equilibrium  of  the  horizontal  bar  (Fig.  8bl : 

fj  ♦  P2  ♦  Fj  -  P  P3  ’  fl  *  °'5P  {i1) 


The  force  in  bar  1  is  now  decreasing,  and  it  will  reach  its 
yield  value  in  compression  while  bar  2  is  still  elastic: 

Stage  JL:  -  -Y  ?2  -  15Y  P  -  12Y  P  *  26Y  (55a) 

kv/Y  =  15  kz/Y  **  11  (55b) 

If  the  load  is  Maintained  at  26Y  with  bars  1  and  3  remaining 
plastic ,  the  solution  consists  of  Bqs.  (55a),  kv/Y  -  IS,  and  z  is 
undetermined  with  any  z  >  0  satisfying  both  relevant  inequalities 
in  (49b) .  Therefore  Bqs.  (55)  represent  the  beginning  of  the 
yield-point  solution. 

Instead  of  eaintaining  P  at  26Y,  let  us  immediately  reduce 
it.  In  moat  problems  a  reversal  of  the  only  load  will  cause  all 
plastic  bars  to  become  elastic.  However,  if  we  superimpose  the 
elastic  solution  (50)  with  a  negative  load  increment  on  Bqs.  (55), 
the  force  in  bar  1  will  immediately  exceed  its  compressive  yield 
strength.  Therefore,  even  though  bar  1  has  just  reached  yield 
under  an  increasing  load,  it  will  resuiin  at  yield  when  we  decrease 
the  load.  The  resulting  solution  is  easily  found  to  be 


Stage  4 :  Pj  -  -Y  P2  -  0.5P  +  2Y  Pj  -  0.5P-Y 

kv  -  2Y  ♦  0.5P  kz  -  11Y 


(56) 


Clearly  no  further  yielding  takes  place  as  P  is  reduced  to 
zero.  When  the  truss  is  fully  unloaded  it  has  a  set  of  residual 
forces  and  permanent  displacements  given  by 

Stage  4L:  P-0  Pj  *  Pj  •  -Y  P2  -  2Y  kv/Y  -  2  kz/Y  -  11  (57) 


The  behavior  of  this  seemingly  simple  structure  becomes  even 
more  complex  when  a  strainhardening  material  is  considered,  under 
an  increasing  load  the  elastic  behavior  and  Stage  1L  will,  of 
course,  be  the  same.  In  Stage  2  bar  1  will  be  in  plastic  tension, 
and  the  stage  will  end  when  bar  3  reaches  its  yield  torce  of  12Y. 

Bar  1  will  unload  and  bar  3  will  be  in  rustic  tension  in  Stage  J. 
Thus  far  the  solution  will  be  only  slightly  different  from  the 
perfectly  plastic  truss  because  of  theeffects  of  strainhardening 
when  bars  1  or  3  are  plastic. 

The  termination  of  Stage  3  will  depend  upon  the  type  of 
hardening.  Por  kinematic  hardening  bar  1  will  yield  in  com¬ 
pression  when  the  change  in  Pj  from  its  value  at  Stage  2L  is  -2Y 
which  occurs  at  a  load  P  -  28. 9Y.  On  the  other  hand,  for  an 
isotropic  hardening  material  bar  1  will  yield  when  its  torce 
reaches  the  negative  of  its  previous  maximum  value  at  P  -  ii.2. 

Pigure  9  shows  the  relation  between  the  load  P  and  the 
deflection 

6  -  v  -  z/2  (5«M 

of  the  point  of  load  application.  The  unloading  portion  shows 
unloading  from  Stage  1L.  Part  of  the  difference  between  the 
curves  is  caused  by  the  fact  that  they  show  unloading  from  different 
values  of  the  load. 

In  the  presence  of  stsinhsrdening  there  is  no  maximum 
allowable  load.  Therefore,  instead  of  reducing  the  load  from 
Stage  3L,  let  us  form  Stage  4  by  continuing  to  increase  it.  bar  1 
will  be  in  plastic  compression,  bar  3  In  plastic  tension,  and 
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bar  2  will  be  still  elastic  with  increasing  force.  The  stage  will 
end  when  F^  reaches  Its  maximum  value  of  20Y.  Table  7  shows  the 
resulting  solutions  for  both  types  of  hardening  at  the  various 
limit  stages,  assuming  r  «  0.1. 

In  Stage  5,  bar  2  will  be  plastic.  If  all  bars  were  plastic, 
the  bar  force  increments  would  be  in  the  same  proportion  as  in 
Stage  1  when  all  bars  were  elastic.  In  particular,  AF1  would  be 
positive  which  would  mean  that  it  would  no  longer  be  in  plastic 
compression.  Therefore,  in  Stage  5  bar  1  will  again  be  elastic, 
while  bars  2  and  3  are  in  plastic  tension.  Stage  5L  occurs  when 
bar  1  reaches  its  current  tensile  yield  stress.  For  the  KH 
material  this  will  be  when  AFj  -  2Y,  whereas  for  an  IH  material 
it  is  not  until  Fj  reaches  the  same  magnitude  tensile  force  it 
had  in  compression  st  Stage  4L.  As  shown  in  Table  7  and  Pigure  10, 
the  load  and  displacement  at  Stage  5L  are  quite  different  for  the 
different  types  of  hardening. 

However,  this  difference  ie  really  qualitative  but  not 
quantitative.  To  see  this,  let  us  continue  to  increase  the  load 
to  a  value  of  SOY.  During  this  Stage  6  all  bare  are  plastic  in 
tension.  Table  7  shows  that  at  the  final  load  all  bar  forces  and 
displacements  have  only  very  small  differences,  and  the  two 
curves  are  virtually  indistinguishable  in  Pig.  10. 

Let  us  review  the  behavior  of  bar  1  of  this  simple 
strain-hardening  truss.  Ae  the  single  load  P  is  monotonically 
increased,  bar  1  is  first  elastic,  then  yielde  in  plastic  tension, 
then  unloads  and  yields  in  plastic  compression,  than  reloads,  and 
finally  yields  for  a  second  time  in  plastic  tension. 
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To  conclude  this  section  we  consider  one  more  loading  program 
in  which  P  is  Increased  until  bar  2  just  reaches  yield  and  is  than 
decreased  to  zero.  Stage  4L  is  the  same  as  before,  but  a  new 
Stage  5*  will  result  from  the  unloading.  Since  fully  elastic 
unloading  would  decrease  f ^  which  is  already  at  compressive 
yield.  Stage  S*  will  find  bar  1  plastic  and  bars  2  and  3  elastic. 

But  the  original  Stage  5  also  had  bar  1  in  plastic  compression. 

In  other  words,  once  the  truss  is  st  Stage  4L,  bar  1  will  continue 
to  yield  in  plaetic  compression  whether  the  load  is  increased  or 
decreased! 

6.  NON -UNIQUENESS 

We  return  to  the  truss  of  Pig.  1,  with  the  following  modi¬ 
fications  the  bare  have  unequal  yield  forces  with  Y^  ■  Yj  •  Y, 

Y2  ■  3Y .  We  consider  e  load -con trolled  program  with  H  -  0  and 
V  Increasing  to  the  yield-point  load.  However,  we  do  not  make 
any  assumptions  of  symmetry.  The  bars  are  all  E/PP. 

Por  convenience  we  shell  repeat  the  defining  equations  in  a 
form  specialized  to  the  particular  loading  program.  The  equili¬ 
brium  equations  are 

fj  ♦  Pj  ♦  /IP2-/JV  <59a) 

Pj  -  Pj  ■  0  159b) 

It  turns  out  that  bar  2  is  always  elastic  end  bars  1  and  3  are 
either  elastic  or  in  plastic  tension.  Therefore,  we  may  combine 
the  klnesMtic  and  constitutive  equations  to  obtain 
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P2  »  kv 

(59c) 

EITHER  Fj  - 

(k/2) (v  ♦  u) 

OR  Fl  •  Y 

{ 59d) 

EITHER  Fj  - 

(k/2) (v  -  u) 

O 

X 

■s 

• 

(59e) 

together  with  inequality  condition*  which  say  be  written 

Fj  i*  P2  <  v  |u|  <  V  (59 f) 

Equation*  (59)  provide  a  total  of  five  equation*  to  determine 
Fi'  P2'  P3'  u*  and  v- 

Stage  1,  of  course,  i*  elastic,  so  we  u*e  the  first  branch 
in  Eqs.  (59d,  e)  to  obtain  the  unique  solution 

u  -  0  2Ft  -  2Fj  -  fj  -  kv  «  (2  -  /I)V  <*01 

as  in  Sec.  2.  However,  with  the  stronger  bar  2,  bars  1  and  3 
yield  first,  so  that  Stage  1L  Is 

u  -  0  fl  *  P3  “  Y  P2  -  kv  -  2Y  V  -  (2  ♦  /I)Y  (61) 
In  Stage  2  bars  1  and  3  are  both  plastic  so  the  second  branch 
is  used  in  (S9d,  e) .  However,  this  Beans  that  (59b)  is  satisfied 
identically.  He  can  still  use  (59a,  c)  to  obtain  the  unique  values 

P1  “  p3  “  Y  P2  -  kv  -  V  -  /?Y  (62a) 

but  we  cannot  determine  the  value  of  the  horisontal  displacement  u, 
although  f$9f)  does  provide  the  bounds 

ju|  <  AV/k  (62b) 

in  particular.  Stage  2b  at  the  yield-point  load  is 

F1  "  P1  "  Y  P2  *  kv  “  3Y  V  -  ( 3  ♦  /I)  Y  (63a) 


Figure  11  shows  the  unique  solution  at  Stage  it  and  several 
solutions  for  Stag*  2'„.  For  example.  Fig.  lib  shows  the  largest 
allowable  value  of  u,  u  »  Y/k.  In  this  solution  bar  3  rotates  a 
a  rigid  body,  bar  2  elongates  elastically,  and  bar  1  elongates 
plastically.  Mot ice  that  any  value  of  u  larger  than  Y/k  would 
require  a  shortening  of  on  3  which  is  not  permissible  when  it 
is  yielding  in  tension.  On  the  other  hand  any  position  of  point 
D  between  the  limiting  one*  in  Figs.  11  a  and  b  requires  a 
lengthening  of  both  bars  1  and  3  which  is  consistent  with  their 
both  being  at  tensile  yield. 

Mathematically,  the  non -uniqueness  was  caused  by  the  fact 
that  two  bars  yialded  simultaneously  and  reduced  one  of  the 
equilibrium  equations  to  an  identity.  In  reality,  of  course, 
infinitesimal  differences  between  the  two  bars  would  make  it 
extremely  unlikely  that  they  would  both  yield  at  exactly  the 
same  instant.  To  examine  this  facet  of  the  problem,  let  us  make 
perturbation  of  the  problem  by  taking  Yj  «  Y  ♦  X  where  X  is 
I  positive  but  otherwise  arbitrary.  Stages  l  and  It.  will  be  the 

same  as  before,  although  only  bar  1  will  yield  at  Stage  1L.  In 
Stage  2  we  now  uee  the  firet  branch  of  (59e)  along  with  the 
second  branch  of  (S9d).  The  resulting  unique  solution  consists 
of  (62a)  tc  -ther  with  the  value 

u  -  4VA 

In  particular,  at  Stage  2l  the  solution  consists  of  Eqs.  (63a) 
plus 


|u|  <  Y/k 


(63b) 


u  -  Y/k 
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Piqure  lib  shows  the  resulting  configuration.  This  solution  is 
completely  independent  of  X,  provided  only  that  X  is  positive. 

In  particular,  X  may  be  arbitrarily  saiall  so  that  there  is  no 
practical  difference  between  the  two  yield  forces. 

However,  it  is  clear  that  if  we  leave  Y^  -  Y  and  set 
Y |  *  Y  ♦  X ,  the  solution  at  Stage  2 L  will  consist  of  Eqs.  163a) 
plus 

u  -  -Y/k  (66) 

as  pictured  in  Fig.  11a.  This  solution,  too,  is  vslid  for 
arbitrarily  small  x.  Therefore,  although  the  non-unique  solution 
of  the  original  problem  can  be  made  unique  by  en  arbitrarily  amall 
perturbation,  two  different  perturbations  give  two  very  different 
unique  solutions. 

Let  us  now  consider  the  effect  of  etrelnhsrdening .  In  doing 
so  we  shall  allow  for  different  rates  of  hardening  in  bars  1  and  3 
by  defining 

V  -  *3'  "  <67) 

Then  Eqs.  (59d,  e)  will  be  replaced  by 

IP  Pj  <  Y  THEN  rx  -  k(u  ♦  v)/2  ELSE  Fj  -  r^lu  ♦  v)  /2  (68a) 

IP  Pj  <  Y  THEN  P3  -  k(v  -  u)/2  ELSE  Fj  •  Tjklv  -  U) /2  (61b) 

Stages  1  and  1l  are  the  same  as  before,  and  the  complete 
solution  for  Stages  2  and  2L  is  easily  obtained.  Unique  values 
are  obtained  for  all  variables,  and  if  the  strainhardenlng  tends 
to  zero  the  bar  forces  and  vertical  displacement  will  tend  to  the 
values  in  Eqs.  (62a)  and  (63s).  Na  shall  focus  our  attention  on 
the  unique  value  of  u  at  Stage  2L  which  Is  given  by 


If  the  strainhardening  in  the  two  bars  is  the  same,  -  r } . 
we  obtain  the  symmetric  eolution  u  -  0  as  pictured  in  Flq.  11c. 

On  the  other  hand.  If  -  0,  and  only  bar  3  hardens,  we  qet 
u  -  -y/k  aa  in  Bq.  (66)  and  Fig.  11a.  Likewise,  it  only  bar  1 
hardens  the  solution  is  again  Eq.  (65),  Fig.  lib.  Clearly 
unequal  non-*ero  hardening  in  the  two  bars  can  produce  any  value 
between  these  two  extremes]  for  example  «  0.1,  Tj  -  0.05  gives 
the  result  u  »  Y/3k  as  shown  in  Fig.  lid. 

7.  CONCLUSIONS 

In  the  preceding  sections  we  have  used  two  very  simple 
trusses  to  illustrate  several  Important  concepts  in  plasticity. 

The  ideas  presented  ere  not  new,  of  course,  nor  are  most  of  the 
applications.  The  trues  in  Fig.  1  has  been  used  many  times  starting 
at  least  as  early  as  1946  [1,  2,  3)*  to  illustrate  plastic  and 
other  Inelastic  behavior.  It  has  also  been  used  to  I lluatrate  some 
singular  features  in  plastic  design  (4) .  The  truss  in  sec.  5 
(Fig.  8}  was  first  introduced  by  Orucker  (5],  and  much  of  the 
present  development  was  taken  from  that  reference. 

Countless  texts  on  plasticity  are  available  (6,  7,  8,  9,  10) 

(to  naam  just  a  few).  The  basic  constitutive  equations  pre¬ 
sented  in  Sees.  1  and  2  can  be  generalised  to  two  and  three 
dimensions  and  to  various  structural  problems.  For  the  particular 

•  Numbers  in  square  brackets  refer  to  references  collected  at 
the  end  of  the  paper. 
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case  of  a  simply-supported  circular  plate  general isations  of  the 
four  stress-strain  curves  shown  in  Pig.  2  have  been  applied  till. 
Figure  12,  taken  f roe  Ref.  [11)  shows  the  relation  between  the 
pressure  and  the  diaplaceeent  of  the  plate  center.  The  results  are 
certainly  qualitatively  sieilar  to  those  for  the  three-bar  truss 
as  shown  in  Pig.  3. 

Shakedown  was  first  introduced  by  Helen  1 1 2 1  .  Many  results 
have  been  presented  by  Symonde  and  by  Neal  (13,  14,  10).  Important 
theoretical  work  has  been  done  by  Koiter  [15,  16).  Recently,  an 
entire  book  by  Gokhfeld  and  Cherniavsky  (17)  has  been  devoted  to 
the  subject. 

The  idea  that  plastic  etreaa  reversal  can  occur  even  under 
a  single  monotonic  loading  has  many  important  implications.  For 
example,  the  somewhat  simpler  theories  known  as  "plastic-deformation” 
theories  which  directly  relate  stress  to  strain  rather  than  relating 
their  rates  is  obviously  inappropriate  wnen  this  phenomenon  occurs. 

Plgure  13  shows  an  eleven-bar  truss  which  does  not  look  par¬ 
ticularly  unusual.  The  top  four  bars  all  have  the  same  cross- 
section  and  are  made  of  an  B/KH  material  with  r  -  0.1.  However 
their  yield  forces  are 

Yj  -  T  Y2  -  Yj  -  12Y  Y4  -  5Y  (70) 

The  bottom  bars  are  all  substantially  stiffer  and  have  a  yield 
stress  high  enough  to  prevent  yielding.  As  the  load  p  is  in¬ 
creased,  bar  1  first  yields  in  tension,  then  bar  4  reaches  its 
yield  force  Sy.  Bar  1  must  now  unload,  and  it  then  yields  in 
compression.  Next,  bar  2  yields  in  tension,  and  bar  1  again 
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unloads,  goes  into  tension,  and  eventually  yields  again  in  tension. 
Table  8  shows  the  bar  forces,  load,  and  displacement  of  the  point 
of  load  application  for  the  following  specific  values: 

Bars  1-4 t  A  “  100  mm2  B  -  80  GPa  fc*  -  8  GPa 

Bars  5-11:  A  -  1000  am2  B  -  240  GPa  (71) 

Y  -  10  MPa 

A  similar  stress  reversal  has  been  observed  in  elastic-plastic 
torsion  of  seme  hollow  bars  (18,  19 1  .  Shaw  (20)  has  nusterically 
solved  the  torsion  of  a  bar  whoae  cross  section  is  a  hollow 
rectangle  with  fillets.  He  has  shown  that  plastic  behavior  starts 
on  the  inner  boundary  at  the  fillet,  and  that  the  stress  vector 
there  is  in  the  direction  of  the  torque.  However,  it  was  first 
pointed  out  in  (18)  and  later  numerically  verified  in  (19)  that 
at  the  limiting  plastic  torque  the  stress  vector  there  must  be  in 
the  opposite  direction  from  the  torque.  An  unfortunate  result 
of  this  plastic  stress  reversal  is  that  the  Nadai  sandhi ll-soapf ilm 
analogy  (21)  is  no  longer  applicable. 

The  phenosenon  of  non-uniqueness  was  encountered,  apparently 
for  the  first  time,  in  the  development  of  a  f inite-element  model 
with  discontinuous  displacements  for  use  in  plane-strain  plasticity 
[22,  23).  Since  this  was  a  primarily  numerical  development,  it 
was  not  clear  whether  the  non -uniqueness  was  inherent  in  the  E/PP 
model,  or  whether  it  was  a  peculiarity  of  the  particular  finlte- 
elsmsnt  model.  The  analysis  presented  in  Sec.  6  (24,  25)  shows 
that  it  is,  indeed,  a  possibility  which  must  be  considered  in  any 
problem  using  the  B/PP  model. 
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Figures  14  and  15  f 24 1  show  some  other  simple  examples.  The 
elastic  solution  for  the  truss  In  Fig.  14  is  symmetric.  It  ends 
when  the  tw"  side  bars  reach  yield.  Further  increase  of  load  can 
be  associated  with  any  non-negative  change  In  the  lengths  of  the 
vertical  bars.  Figure  14  shows  the  two  extrstae  positions  where  one 
of  the  plastic  bars  does  not  change  Its  length.  Any  intermediate 
solution  can  be  obtained  by  a  rigid-body  rotation  of  the  upper 
triangle  about  the  center-point  of  the  truss. 

As  the  load  on  the  frame  in  Pig.  15  is  increased,  hinges 
will  form  first  at  C,  then  at  B,  end  then  simultaneously  at  points 
A  and  d,  but  the  frame  will  not  collapse  until  the  final  hinges 
form  at  E.  The  solution  is  unique  during  the  elastic  and  firet 
two  partly-plastic  stages,  but  once  the  hinges  form  at  A  and  D, 
the  rotations  at  these  new  hinges  are  controlled  only  by  Inequa¬ 
lities,  both  heinq  required  to  exhibit  tensile  strain  at  the  top 
of  the  hinge.  Figure  IS  shows  the  two  extras*  positions.  Any 
intermediate  solution  obtained  by  a  rigid-body  vertical  motion 
of  the  central  part  of  the  frame  is  also  possible. 

The  possibility  of  part  of  the  solution  balng  non-unique 
has  Important  implications  in  finite-element  programs.  Even  If 
the  engineer  is  concerned  only  with  forces  which  ars  unique,  a 
non-unique  soluticn  corresponds  to  a  singular  sat  of  aquations, 
and  hence  to  a  sinqular  stiffness  matrix.  Thus  it  may  be 
impossible  to  continue  the  solution  past  the  point  whets  the 
non -uniqueness  occurs,  evert  though  the  load  is  still  below  the 
yield-point  load.  Mimerlcel  rcund-off  error  may  sask  the 
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singularity  of  the  stiffness  matrix,  but  the  resulting  problem 
would  be  highly  ill-conditioned  and  might  lead  to  misleading 
results.  This  problem  is  certainly  worthy  of  further  investigation. 
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undeformed  position 
uniquely  determined  deformed  position 
possible  non- unique  deformations 
yield  hinges  prior  to  collapse 
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